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ABSTRACT 

We calculate the nucleons' electromagnetic polarizabilities in heavy baryon chiral per- 
turbation theory including all terms to order 0{q' i ). The chiral prediction of the electric 
polarizabilities for the neutron and the proton are in good agreement with the data. In 
the case of the magnetic polarizabilities the big positive contribution from the A(1232) 
resonance is largely cancelled by a non-analytic loop contribution of the lnM x type. 
This novel effect helps to understand the rather small empirical value of the nucleons' 
magnetic polarizability. 



BUTP-93/22 August 1993 

CRN 93-39 

* Work supported in part by Deutsche Forschungsgemeinschaft and by Schweizerischer 
Nationalfonds. 

' Heisenberg Fellow. Address after September 1 st , 1993: CRN, Physique Theorique, 
BP 20, F-67037 Strasbourg Cedex 2, France. 







1. The electric (a) and magnetic (j3) polarizabilities of the proton and the neutron 
encode information about the structure of these particles in the non-perturbative regime 
of QCD. Over the last few years, very precise experiments supplemented with dispersion 
sum rules have established rather tight bounds on these quantities. The outstanding 
features of these results are that the proton and the neutron both behave essentially 
as (induced) electric dipoles (a p ~ d n ^> (3 P ~ (3 n ) and that the sums of electric and 
magnetic polarizabilities are equal within ten percent, (a + f3) p ~ (a + f3) n . The first of 
these observations appears rather stunning since the very strong magnetic NA transition 
gives rise to a large positive contribution to the magnetic polarizabilities. Therefore, it 
has been argued by many that a theoretical explanation of these fundamental quantities 
has to involve explicit A(1232) degrees of freedom. This theme will be picked up below. 
On the theoretical side, chiral perturbation theory (CHPT) offers the best venue of 
understanding these structure constants. Indeed, in refs.[l,2] it was shown that the 
nucleon electromagnetic polarizabilities to leading order in the chiral expansion are 
pure one-loop effects. Consequently, at this order they can be calculated without any 
unknown parameters. This is similar to the prediction for the charge radius of the 
neutral kaon or other mesonic processes like 77 — > 7r°7r°, Kl — > 7r°77 and Ks — + 77. In 
refs.[l,2], the nucleon was treated as a fully relativistic Dirac field. The fact that the 
nucleon mass does not vanish in the chiral limit poses some problems to the chiral power 
counting in the effective meson-nucleon Lagrangian. In particular, the one loop result 
leads to a string of terms in increasing powers of the pion mass. From these, only the 
first one is not affected by contributions from higher loop orders not yet calculated. In 
fact, the calculation of refs.[l,2] was redone in ref.[3] in the framework of heavy baryon 
CHPT (HBCHPT), which allows for a consistent chiral power counting, i.e. loops are 
suppressed by powers of q 2 with q a genuine small four-momentum (this is not the case 
if one treats the nucleon fully relativistically) . In the framework of HBCHPT, only 
the leading singular term proportional to 1/M n in the electromagnetic polarizabilities 
survives and one finds a surprisingly good description of the data (see below). However, 
as mandated by the decoupling theorem [4], in this calculation the intermediate states 
are only nucleons, i.e. no A's appear. As already mentioned, the absence of the low-lying 
spin-3/2 resonances in this approach has prompted many critical remarks concerning the 
results to leading order 0(q 3 ). Jenkins and Manohar [5] have argued that it is mandatory 
to include the spin-3/2 decuplet in the effective theory to cancel large SU(3) loop effects. 
In this spirit, Butler and Savage [6] have performed a calculation of the polarizabilities in 
SU(3) 1 ^ with explicit decuplet degrees of freedom. However, the whole approach suffers 

^ Notice that the first SU(3) calculation of these observables was performed in ref.[7]. 
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from two shortcomings. First, including the decuplet destroys the consistent power 
counting of HBCHPT. This can be traced back to the non-vanishing scale related to 
the finite decup let-octet mass splitting which does not vanish in the chiral limit [8]. 2 ) 
Furthermore, the calculations performed in refs.[5,6] do not include all terms in a given 
order in CHPT but only some non-analytic loop contributions and some analytic terms. 
Clearly, this does not conform to the philosophy underlying CHPT, namely to include all 
terms at a given order. The A(1232) was also hailed as a deus ex machina in explaining 
the data of the nucleon electromagnetic polarizabilities by L'vov [9]. In this paper, we 
will present a complete CHPT calculation of these fundamental quantities including all 
effects up to (and including) order q 4 instead of arguing about this or that being the 
dominant effect. 

2. Consider the spin-averaged Compton amplitude for scattering off nucleons at low 
energies. The energy expansion of the spin-averaged Compton amplitude in the nucleon 
rest frame takes the form 

e 2 Z 2 - 
T("fN -> 7 JV) = + au'ue' ■ e + (3 (e' x k') ■ (e x k) + 0(u 4 ) (1) 

with k, e) and (u/, k! ', e') the frequencies, momenta and polarization vectors of the in- 
coming and outgoing photon, respectively. The first term, which is energy-independent, 
is nothing but the Thomson scattering amplitude as mandated by gauge invariance. 
It constitutes a low-energy theorem as the photon energy tends to zero. At next-to- 
leading order in the energy expansion, the photon probes the non-trivial structure of 
the spin-1/2 particle (here, the nucleon) it scatters off. This information is encoded 
in two structure constants, the so-called electric (a) and magnetic (/?) polarizabilities. 
These have been determined rather accurately over the last years. For the proton, if 
one combines the Illinois, Mainz and Saskatoon measurements [10] with the dispersion 
sum rule, (a + f3) v = 14.2 ± 0.3 [11], one has 

a p = 10.4 ±0.6, P = 3.8^0.6 (2) 

all in units of 10 -4 fm 3 which we will use throughout and do not exhibit explicitely any 
further. Notice that we have added the statistical and systematic errors in quadrature. 

2 ) This might be circumvented if one considers the NA mass splitting as a small param- 
eter. However, in three-color QCD this quantity can not be considered arbitrarily small 
in contrast to the light quark masses or external momenta of the fundamental fields in the 
effective field theory. 
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Similarly, the dispersion sum rule (a + (3) n = 15.8 ± 0.5 [12] 3 ) together with the recent 
Oak Ridge and Mainz [13] measurements lead to 

a n = 12.3 ± 1.3, p n = 3.5 =f 1.3 (3) 

The salient features of these experimental results have already been mentioned, namely 
that both the proton and the neutron essentially behave as (induced) electric dipoles and 
that their respective sums of electric and magnetic polarizabilities are almost the same. 
Let us stress again that naively one expects a large contribution from the A(1232) to 
the magnetic polarizabilities (of the order of 10) and it thus appears difficult to explain 
the relatively small values of j3 Pjn . 

3. To calculate the polarizabilities, we make use of HBCHPT. We work in flavor SU(2) 
and in the isospin limit m u = ma- The relevant degrees of freedom of the effective 
Lagrangian are the pions and the nucleon. The Goldstone fields are collected in the 
matrix U(x) = exp[zr •jf(x)/F 7T ] = u 2 (x) and the nucleons are considered as very heavy, 
i.e. non-relativistically [14]. The effective Lagrangian to order 0(g 4 ), where q denotes a 
genuine small momentum or a meson (quark) mass, reads (we only exhibit those terms 
which are actually needed in the calculation) 

r - jl r( 2 ) jl X 4) 
ci 1} N = H(ivD + g A S-u)H 

4 2 ^^{-^^ + ^^) 2 + Cl Tr X+ + ( C2 -^.^.. 

+ C3 u ■ u - ^{S ■ D, v ■ u} - S»] ((1 + «„)/+ + Tr /+ ) }ff (4) 

4* = ~ Wn)Hf+J+ H + jS0 n HHTrf+ 



+ ^(da n + Sp n - da p - 6P P )Hf+J$ v Hv' t vx 
- ^(5a n + 5p n )HHTr(f+J^)v»v x 



with 

(5) 



3 ) Notice that the error on this number is presumably underestimated since one has to 
use deuteron data to extract the photon-neutron cross section. 
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where H denotes the heavy nucleon field of charge Z = (1 + ts)/2 and anomalous 
magnetic moment n = (k s + tsk v )/2, the covariant spin-operator subject to the 
constraint v ■ S = 0, the covariant derivative acting on the pions and we adhere to 
the notations of ref.[3]. The superscripts (1,2,4) denote the chiral power. The lowest 
order effective Lagrangian is of order O(q). The one loop contribution is suppressed with 
respect to the tree level by q 2 thus contributing at 0(q 3 ). In addition, there are one loop 
diagrams with exactly one insertion from C^' N . These are of order q 4 . Finally, there are 
contact terms of order q 2 and q 4 with coefficients not fixed by chiral symmetry. For the 
case at hand, does not contribute. Notice that some coefficients in C^ N related to 
the *yyNN and 'jirNN vertices are fixed from the relativistic theory by the low-energy 
theorems for Compton scattering and neutral pion photoproduction, respectively. This 
is discussed in some detail in ref. [3] . The unknown coefficients we have to determine are 
ci, C2 and C3 as well as the four low-energy constants 5a p , 5a n , d(3 p and 5j3 n from C^ N . 

(2) 

We have not exhibited the standard meson Lagrangian L\i . 

To calculate all terms up-to-and-including order g 4 , we have to consider all one 
loop graphs with insertions from and those with exactly one insertion from 
While the former scale as o 3 , the latter constitute the new contributions of order q 4 . 
Furthermore, there are the tree diagrams related to C^ N which are also new. It is worth 
to stress that in the one loop diagrams involving C\. ^ the anomalous magnetic moment 
of the nucleon appears. Finally, it is mandatory to expand the leading order O(q) 
effective vertices as well as the nucleon propagator in momenta to include all relativistic 
corrections of order 1/m. The details of these calculations can be found in ref. [15]. For 
further convenience, we introduce the constant C, 

c = s^?f = 4 - 36 ' 10 " 4fm2 <6) 

with e 2 /4n = 1/137.036 and F n ~ 93 MeV the weak pion decay constant. The result- 
ing expressions for the electric polarizabilities of the proton and the neutron can most 
compactly be written as 



_ 5Cg 2 A C 
ai ~ 4M^ + n 



(^-c 2 )ln^ + i(MA-6c 2+C +) 
v m ' A 4 2m ' 



+ Sa r i (X), i = p,n (7) 



with 

x p = 9 , x n = 3 , y p = 71, y n = 39. (7a) 

Here, A is the scale of dimensional regularization. The physical oti are of course indepen- 
dent of this scale since the renormalized counter terms dd^(X) cancel the logarithmic 
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A-dependence of the loop contribution. The first term on the r.h.s. of eq.(7) is, of 
course, the result at order q s [3]. Similarly, one finds for the corresponding magnetic 
polarizabilities 



8M OT 



TV 



(5£M_ C2)ln ^ + I ( M + 2c2 _ c+) 

m ' A 4 m ' 



37 



+ 6k s 



Vn 



13 



with 

We have introduced the coefficient c + related to ci, c 2 and C3, 



- 6k, 



p, n (8) 



(8a) 



,+ _ 



9a 



= -8ci + 4c 2 + 4c 3 - ^ 

2m 



(9) 



The results shown in eqs.(7,8) have the following structure. Besides the 1/M n term, 
which is discussed in detail in refs.[2,3], the O(o 4 ) contributions from the loops have a 
In M n and a constant piece ~ M®. As a check one can recover the coefficient of the In M n 
term form the relativistic calculation [1,2] if one sets the new low energy constants Cj 
and k SjV = 0. In that case only the 1/m corrections of the relativistic Dirac formulation 
are treated and one necessarily reproduces the corresponding non-analytic (logarithmic) 
term of this approach. The term proportional to c 2 lnM^ in eqs.(7,8) represents the 
effect of (pion) loops with intermediate A(1232) states [6] consistently truncated at 
order g 4 . Notice that from now on we will omit the superscript V on 5a^ and 8(31 
appearing in eqs.(7,8). Finally, from eqs.(7) and (8) one can form two quantities which 
are independent of the low-energy constants c 2 and c + and which vanish to leading order 



£+ = (a + /3) p -(a + f3)i 

= Cg\L,„ M„ 
tv m 

£- = («- 0) p -(a-P) 



6(2 + Ks )ln^ + (7 + 3« s ) 
A 



+ Sa p (X) + Sp p (X) - Sa n (X) - Sp n (X) 



(10) 



tv m 



2 v 



l-3 Ks (l + 21n— -) 
A 



+ 5a p (X) - Sp p (X) - 5a n (X) + 8(3 n {X) 



Empirically, one has E + = —1.6 and E~ = —2.2. From these numbers one immediately 
realizes the importance of isospin breaking in the counter terms from C^ N . Setting their 
renormalized values equal for the proton and the neutron and choosing A = 1.232 GeV, 
one finds E + = —8.3 which deviates considerably from the empirical value. Similarly, 
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the prediction E = —0.1 is not in good agreement with the data. We will come back 
to these quantities later on. 

4. We now have to determine the values of the various low-energy constants appearing 
in eqs.(7) and (8). The constants 01,2,3 already appeared in the study of the chiral 
corrections to the S-wave pion-nucleon scattering lengths [16]. From that paper, we 
can deduce the relation 

C ~ 32nM^(l-M 7T /m)- {U) 

where a + is the isospin-even S-wave scattering length, a+ = —0.83 ± 0.38 • 10 _2 /M 7r , 
gA the axial-vector coupling constant and M n = 139.57 MeV the (charged) pion mass. 
Inserting the value of c + as given from eq.(10) into eqs.(7,8), one finds that the con- 
tributions proportional to c + to the electromagnetic polarizabilities are negligible since 
they are less then 0.12 in magnitude. To determine C2, we make use of the principle of 
resonance saturation [17]. It states that to a high degree of accuracy the low-energy 
constants can be calculated from resonance exchanges by integrating out the heavy res- 
onance fields from an effective Lagrangian of the pions chirally coupled to the various 
resonances. In the meson sector, this has been shown to work very well. We extend this 
method to the baryon sector since it is essentially the only method of estimating the 
unknown coefficients. For the constant C2 we find 

C2 = 0( 9 " A ; =0-59fm (12) 
2(?tia — Tn) 

where we have treated the A non-relativistically (isobar model) and neglected a small 
contribution from the Roper resonance. We use wia = 1232 MeV and m = 938.27 MeV. 
A detailed discussion on the admissible range for c<i including also its dependence on the 
off-shell parameter Z if the A is treated relativistically (as a Rarita-Schwinger spinor) 
can be found in ref.[15]. Here, we will use the value given in eq.(12) at the scale A = ira 
according to the resonance saturation principle. 

The A(1232) enters prominently in the determination of the four low-energy con- 
stants from C^pf. Therefore, we will determine these coefficients at the scale A = m&. 
In particular, one gets a sizeable contribution to the magnetic polarizabilties due to 
the strong iVA Ml transition. A crude estimate of this has been given in ref.[18] by 
integrating the Ml part of the total photoproduction cross section for single pion pho- 
toproduction over the resonance region, 

S^(m A ) = ^-f^a M \u) = 7.0 (13) 
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However, there is a large uncertainty this number. If one simply uses the Born dia- 
grams with an intermediate point-like A, one has 5(3^ (m A ) = 10 ... 12 [18] where the 
uncertainty stems from the uncertainties in the 7./VA coupling strength. It has been 
argued in ref.[19] that an off-shell effect could reduce this number by almost an order of 
magnitude, however, this is outside the accuracy of the chiral expansion performed here. 
Furthermore, from isospin arguments one can conclude that 5(3^ (ttta) = 5/3^ (toa) — 7. 
In addition to this resonance contribution, one has also an effect from the (charged) 
kaon loops [7]. Since we are working in SU(2), the kaons and etas are frozen out and 
effectively give some finite contact terms. One has for the electric polarizabilities 

K 5C 2D 2 /3 + 2F 2 

Sa p ^ = x ~Mk = 

K (14) 

for D = 0.8 and F = 0.5. The corresponding numbers for the kaon contributions to the 
magnetic polarizabilities are a factor 10 smaller [7]. With this, we have determined all 
the unknown coefficients and are at the position of evaluating eqs.(7,8). A more detailed 
discussion of these issues can be found in ref.[15]. 

5. We now present our numerical results. We always use the Goldberger-Treiman 
relation to express qaIF^ as g^N /m , with g n N = 13.40 the strong pion-nucleon coupling 
constant. We find with k s = —0.12 

a p = 13.6 - 8.3 + 3.2 + 2.0 = 10.5 
a n = 13.6 - 1.6 + 1.2 + 0.2 = 13.4 

(15) 

P P = 1.4 - 7.9 + 2.8 + 7.2 = 3.5 
p n = 1.4 - 2.0 + 1.4 + 7.0 = 7.8 

where we have separated the q 3 contribution (first term) from the q 4 ones which are of 
the lnM,,-, M° and counter term type. Several remarks on these results are in order. For 
the electric polarizabilities, we find small corrections which bring the chiral prediction 
in agreement with the data, in particular the difference in the proton to neutron electric 
polarizabilities. For the magnetic polarizabilities, the most important result is that the 
large contribution from the A(1232) which enters through the counter terms is mostly 
cancelled by the large coefficient of the In term in case of the proton. For the neutron, 
this term is smaller and therefore the resulting chiral prediction is too large. Clearly, 
to overcome this at order g 4 , one would need some substantial isospin breaking in the 
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photon-nucleon-Delta coupling constants (for which at present we have no indication). 
Consequently, the quantity E + defined in eq.(10) comes out too big, E + = —7.2. Also, 
we find E~ = 1.5 which has the wrong sign due to the to large magnetic polarizability 
of the neutron. Is it evident at this stage that a better understanding of the counter 
terms from C^ N is necessary to further tighten the chiral prediction for the magnetic 
polarizabilities of the proton and the neutron. 

6. To summarize, we have used heavy baryon chiral perturbation theory to calculate all 
corrections up-to-and-including order q 4 to the electric and magnetic polarizabilities of 
the proton and the neutron. To estimate the strength of the various contact terms, we 
have made to some extent use of the principle of resonance saturation which is known 
to work very accurately in the meson sector. The main results of this investigation are: 

• The corrections to order q 4 bring the chiral prediction of the proton and neutron 
electric polarizabilities in agreement with the data, d p = 10.5 • 10~ 4 fm 3 and a n = 
13.4-10 -4 fm 3 . The corrections to the leading order (q 3 ) result = 13.6-10" 4 
fm 3 are small and one thus expects no drastic changes to these results from the next 
order, 0(q 5 ). 

• The situation is different for the magnetic polarizabilities. There is undoubtedly a 
large positive contribution from the A(1232) resonance. Its actual value, however, 
is quite uncertain. For the proton, one finds a substantial cancellation between the 
positive A and large negative lnM^ term bringing the chiral prediction in good 
agreement with the data. For the neutron, the lnM x term has a considerably 
smaller coefficient and thus the resulting value for j3 n is approximately a factor two 
too large. 

A more detailed account of these topics, in particular a discussion of the allowed ranges 
for the various counterterms (which was not addressed in detail here) will be presented 
in ref.[15]. Finally, it is worth to stress again that the calculation presented here contains 
all effects up-to-and-including order q 4 in the chiral expansion of the nucleons' electro- 
magnetic polarizabilities and that the resulting numbers reflect most of the (although 
not all) empirical facts about these fundamental quantities. 

We are grateful to Alan Nathan, Nimai Mukhopadhyay and Ron Workman for instruc- 
tive comments. 
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